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1. Introduction 



Let Q be a bounded domain in R 3 and Qt = x (0, T] for any fixed T > 0. Let E 
and H be the electric and magnetic fields, respectively, in Q (here and thereafter a 
bold letter represents a vector in R 3 ). Let a be the electric conductivity in the field, 
which is assumed to be a function of x and |E|. Consider the following Maxwell's 



equations (see Landau-Lifschitz |fL3 



eE t + a(x, |E|)E = V x H + F, (x,t) G Q T , (1.1) 

H, + VxE = 0, (x, t) G Qt, (1.2) 

VxE = 0, (x, t) G <9J1 x (0, T], (1.3) 

E(x, 0) = E (x), H(x, 0) = H (x), x e tt, (1.4) 

where e is the dielectric parameter and other physical parameters are normalized. 

In some applications ([§], [TIH), the electric conductivity, cr, strongly depends on 
the electric field |E|, hence the electric current density. Particularly, the electric 
conductivity may act like a switch-like function in some electromagnetic fields. On 
the other hand, for many types of micron devices and other industrial problems (such 
as microwave heating ]TI, |8|, [H[]) etc.) the experiment shows that the displacement 



current, eE t , is often negligible since it is small in comparison of the eddy current, 
J = crE. This motivates us to investigate the nonlinear problem (1.1)-(1.4) and 
the singular limit problem as e — > 0. It is shown that there exists a unique global 
solution to (1.1)-(1.4). Moreover, the limit of the solution converges to the solution 
of the quasi- stationary system (i.e., the system (1.1)-(1.4) with e = in (1.1)). This 
limit solution provides new existence result for the quasi- stationary system. Indeed, 
when e = 0, the system (1.1) becomes 

a(x, |E|)E = V x H + F, (1.5) 

Thus, one can solve Eq. (1.5) for |E| in terms of |V x H| and known data, 

|E| =g{x,\Vx H|), 
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where g(x, s) is the inverse function of a(x, s)s. 
It follows from (1.2) that H satisfies 



H t + Vx[p(i,|VxH|)VxH] = 0, (1.6) 

where 

p(x,\VxH\\) 



a(x,\E\) a(x,g(x,\V xH\)) 
represents the electric resistivity in the field. 

The research on Maxwell's equations is of great interesting because of the im- 
portant applications in plasma physics, semiconductor-superconductor modeling and 
other industrial problems (0, || |13|, [17] etc.). The study on the system (1.1)-(1.4) 
as well as the quasi-stationary form (1.6) received considerable attention recently. In 
[F2], the authors established the well-posedness for a quasi-stationary system, where 



a constitutive relation between the magnetic field H and the magnetic induction B is 



assumed to be nonlinear. In [flql , the author studied the regularity of weak solution to 
a linear system of (1.6) with minimal requirement on coefficients. There is a special 
interest when 

p(|V x H|) = |V x H| p ~ 2 , p>2. 

On one hand, if H is restricted in one direction (scalar field) then the evolution system 
(1.6) becomes the p-Laplacian which has been studied extensively (see f| and the 
references therein). On the other hand, in a recent work jnj (also see |1], |I| for the 
scalar case), it is shown that the limit of the solution to (1.6) as p — > oo is the unique 
solution to Bean's critical-state model in the superconductivity theory (@). Thus, 
for large p the system (1.6) provides a good approximation to Bean's model. More 
recently, the author of |I| studied the similar problem to this paper in a domain 



with a bounded complement in R 3 . The conditions on a in |10| is quite different 



from ours here. Like many nonlinear problems, the major difficulty is how to pass 
the weak limit of an approximate solution for a nonlinear function a(x, s). This is 
done by employing a monotonicity argument (H). The monotonicity of a(x, s) in s 
is essential in the proof. 
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In §2, we use the finite element method to establish the well-posedness of the 
system (1.1)-(1.4) for fixed e > 0. In §3, we show that the singular limit of the 
solution to (1.1)-(1.4) has a unique limit. Moreover, the limit solution solves the quasi- 
stationary Maxwell's equations. Some examples are also discussed in this section. 

2. Existence and Uniqueness for fixed e > 

Introduce some standard spaces (see 0, 0). 

H(div,fl) = {U G L 2 (£lf,divXJ G L 2 (tt)}; 
H(curl, fl) = {Ue L 2 (ttf, curlXJ G L 2 (ft) 3 } ; 
H(div0, fl) = {Ue H(div, Q) : divU = in ft}, 
H (curl,fl) = {U G H(curl,fl) : N x U = on d£l}, 

where N is the exterior unit normal on dfl. 

Note that the trace of a function in H(curl, Q) is well defined (see 0] for example). 

We shall assume the following conditions on a(x,s) and data Eo(x),Ho(x) and 
F(x,t). 

H(2.1): Let a(x, s) be measurable in Q x [0, oo) and monotone increasing in s. More- 
over, 

/ a(x, \fs)dx > a s p+2 — a 1; if s is sufficiently large, 
Jo 

< a(x, s) < b (l + s p ), s G [0, oo), for some p > 0, 
where the constants a > 0, a\ > and &o > 0. 

H(2.2): Assume that H G H(curl, O) f] H(div0, O), F G -£^(0, T; H 1 ^)). 

Definition 2.1: A pair of vector fields (E(x,t),H(x,t)) is said to be a weak solution 
of the problem if 

E G L 2 (0,T;i/ (aiW,fi))fl^ p+2 (0,T; fi),H G L 2 (0, T; H(div0, SI)) f]H(0, T, L 2 (Sl)) 
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which satisfy the following integral identities: 

f f [-eE-& t + a(x,\E\)E-&]dxdt 

= 11 H • (V x &)dxdt + e [ E ■ <&{x,0)dx, (2.1) 
J Jq t Jn 

[ [ [-H • + E • (V x *)] dxdt = [ [H (x) • *(x, 0)] dx (2.2) 
J JQ T Jn 

for all test functions $ G H l (0,T; H (curl,tt)),V G H\0,T; H (curl, SI)) with 
&(x,T) = *(x,T) = 0. 

First of all, we derive some energy estimates. A special attention is paid on how 
various constants depend on e since we will study the singular limit problem in section 
3. 



Lemma 2.1: Under the assumptions H(2.1)-H(2.2) there exist constants C\,C2 and 
C 3 such that 



sup 

[o,T] Jn 

Jn 



e\E\ 



IHI 



| E |p+2' 



dx + 



e|E t | 



IHJ 2 



io Jn 

2 TJ 2 V7 TT |2 J„ i r< i I NT7|2 IT 2 



|E | 2 + |H | 2 + |V x H | 2 dx + C 2 



dxdt 

F\ 2 + \F t \ 2 ]dxdt + C 3 , 



where C\,C2 and C3 depend only on known data. 

Proof: Note that for any vector fields A, B G H(curl, Q) with either A or B in 
H (curl,fl), the following identity holds: 



J A • (V x B)dx = B • (V x A)dx. 



Taking inner product to the system (1.1) and (1.2) by E and H, respectively, we add 
up the resulting equations to obtain 



In 

<of 

Jn 



sup 

[0,T] Jn 



£|E| 2 + |H| 2 " 


dx + j 


[ a(x, \E\)\E\ 2 dxdt 






JQt 



I Eq 1 2 + Hq I 



dx+ f f [\E-F\]dxdt 
J Jq t 



(2.3) 



where the constant C depends only on known data, but not on e. 
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We first assume that a(x, s) is differentiable with respect to s. Then we formally 
differentiate Eq.(l.l) and Eq.(1.2) with respect to t to obtain 



sE tt + a(x, \E\)E t + a s (x, |E|)(|E|) t E = V x H t + F t , 
H tt + VxE ( = 0. 

It is clear that 

/ / (V x Et) • Utdxdt =11 (V x H t ) • E t dxdt. 

J JQ T J JQ T 

We take the inner product by E t for the first equation and by H t for the second 
equation and add up the resulting equations to obtain: 



/ 

Jn 



sup 

[o,t] 

<C, 



e\E t \ 2 + |Ht| 2 dxdt + 



a\E t \ 2 + (7 s (|E|) t E-E t 



dxdt 



where C depends only on known data. 
Note that a s > 0, we see that 



/ / <r a (|E|) t E • E t dxdt 



d IEI 2 
a s (\E\) t - l -^-dxdt 



I Iq t 



dt 2 
\E\\\E\ A 2 dxdt > 0. 



It follows that 



sup 

[o,T] Jn 



e\E t \ 2 + I H* 



dxdt + / / [(r|Ei| 2 dxdt < C. 



Since the above estimate does not depend on the differentiability of a with respect to 
s, therefore the above estimate holds as long as a is monotone increasing with respect 
to s. 
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Now we take the inner product by E t to (1.1) and by H t to (1.2) and then add 
up the resulting equations to obtain 



< 



1 , r\E(x,T)\* 

dxdt + - / / a(x,^s)dsdx 

2 Jn Jo 

'IEoMI 2 



/ / E t -Fdxdt + C ( ( °* a(x,^)dsdx+ f |VxH | 2 dx + C. 
J JQ T Jn Jo Jn 



Now 



E t ■ Fdxdt 

[ [E(x,T) -F(x,T) -E (x) -F(x,0)]dx- f [ E F, 
Jn J Jq t 



dxdt 



< [ f^E(:r,T)r 2 + ±V t |^ 
Jn L 4 a Q 



dxdt. 



On the other hand, by the assumption H(2.1) we may assume that the growth con- 
dition of a(x, s) on s holds for all s > M , i.e., 



/ a(x, y/s)dx > a s p+2 — cii, if s > M , 
Jo 



where M is a fixed constant. 
It follows that 

-|E(*,T)p 



/ / a(x, s)dsdx 

Jn Jo 

>a [ \E(x,T)\ p+2 dx-C. (2.4) 

Jnft{x:\E(x,T)\>M } 

Combining (2.3)-(2.4) yields 

sup f [e\E\ 2 + \H\ 2 ]dx + sup f \E(x, t)\ p+2 dx + f [ [e\E t \ 2 + \U t \ 2 ]dxdt 

<t<rJn o<t<rJn Jo Jn 

[\F\ 2 + \F t \ 2 ]dxdt + [ [|E | 2 + |H | 2 + |V x H | 2 ]cfe + C. 
Jn 

Q.E.D. 



0<t<T 

< c 
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Theorem 2.2: Under the assumptions H(2.1)-H(2.2) the problem (l.l)-(l.^.) has a 
unique weak solution. Moreover, 

curiE G L 2 ((5r),E t G L 2 (Q T ) 

and 

a t e L 2 (Q r ),V x H g lSt(Q t ). 



Proof: The proof is based on the finite element method (see [14[] for parabolic equa- 
tions). The monotonicity of a(x, s) on s plays an important role. We shall first deal 
with the case where a(x, s) is continuous on s. For convenience, we rewrite the system 
(1.1)-(1.4) to the following form: 

sW tt + a(x, |W t |)W t = V x [H - V x W] + F, (x,t) G Q T , (2.5) 
N x (W ( ) = 0, x G dn, < t < T, (2.6) 

W(x,0) = 0,W t (x,0) = E (a;), x£fi, (2.7) 



where W is defined as follows: 



W{x,t) = /"*E(x,T)dr 
Jo 





It is clear that if W is a solution of the system (2.5)-(2.7) then a pair of functions 
defined by 

E(x, t) = W t (x, t), H(z, t) = H (x) - V x W(ar, t) 

will be a weak solution of (1.1)-(1.4). Let {e^} = {e^ , e^ 2 , ej^} be a smooth basis 
of H (curl,Q) and orthonormal in L 2 (Q) 3 , i.e. 

<C e^ j e j ^> 5{j , 

where <5y = 1 if i = j and ^ = if i ^ j. 

Now we expand the known data as follows: 

oo 

H o(^) = dia 9[*k o e fc ], 
fe=i 



S 



E o(aO = J2diag[b k oej, 
fc=i 

oo 

F(x, t) = 2 dm5f[g fc (t) o e fc ] 



k=l 

where a^, b k and are 3x1 matrices, the symbol o is the matrix product and diag\\ 
represents the diagonal vector of a matrix. 
Let 

n 

W n (x, t) = J2 diag[cW o e fe ], 
k=i 

where c^\t) is a 3 x 1 vector which is determined by the following ordinary differential 
system: 

e|U fe) + \ W ^j/n ) = e fc ) + ifc(f), (2.8) 

C?>)(0) = 0, (2.9) 
|cf(0)=b fc , (2.10) 

where 

WW = dia0[c?>oe fc ], 

4(W«,e fc ) = J Q dia 9iN x WW] o [V x e k ]}dx, 

B k (t) = { diag{(V x H + F) o e k }dx, k = 1, 2, • • • n. 
Jn 

Now we define the approximate solution (E n ,H n ) as follows: 

E n (x,t) =Wra(x,t),H n (x,t) =H 0n - V xW n (x,t), 

where 

n 

H on(^) = diag[a k o e fe ]. 
k=i 

Equivalently, then (E„,H n ) satisfies the following system in the weak sense: 

eE nt + a(x, |E n |)E n = V x H n , (x,t) e Q T , (2.11) 
H nt + VxE n = 0, (x,t)eQ T . (2.12) 



Similar to Lemma 2.1, one can easily derive the following energy estimates : 



sup / 

[0,T] 

Jn 



e\E n \ 2 + |H„| 2 + |E n | p+2 l dx + 



^|E n J 2 + |H 



IE. 



On | 



+ |H J 2 + IV x H 



On I 



o Jn ' 

dx + C 2 j [ [|F r 
J JQt 



nt I 



dxdt 



+ \F nt \ 2 ]dxdt + C 3 , 



where Ci, C2 and C3 are independent of n and e. 

By the weak compactness property, we can extract a subsequence (still denoted 
by (E n ,H„)) such that 

E n -> E, E ni -> E t , H ni -> H t , weakly in L 2 (Q T ), 
H n -> H, weakly in L 2 (0, T; V^ 1 '^ (fi)), 
H n — > H, a.e. in Qr- 



Moreover, 



E r 



E 



weakly in L P+2 (Q T ). 



Next we claim that the sequence E n converges to E strongly in L 2 (Qt) 3 - To prove 
the claim we only need to show that {E n } is a Cauchy sequence in L 2 (Qt) 3 . Let 

E* n {x, t) = E n (x, t) - E m {x, t),H* n (x, t) = H n {x, t) - U m (x, t). 

By energy estimates, we see 

sup /[|e;i 2 + |h;i 2 ]^ + 

f f {[cr(x, |E n |)E„ — a(x, |E m |)E m ] • [E„ — E m ]}dxdt 

J JQt 



<c I 

Jn 



Eon — E 



0m| 



+ IH 



„ - H 0m \ 2 ]dx + C f f [|F n - F m \ 2 }dxdt, 



where C is a constant independent of n and m. 

Note that a(x, s) is monotonic increasing in s, then 



[a(x, |E n |)E n u{x, |E m |)E r 



E„. — E r 



> 



a(x, \E n \) - a(x, \E r , 



E|2 IIP 



> 0. 
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It follows that 



sup / [|e;i 2 + |h;i 2 ^ 



'dxdt. 



This implies that E n , H n are Cauchy sequences since both Eo n , H 0r! , and F n are 
Cauchy sequences in L 2 (Qt) 3 - Hence, 

E n ,H„ -> E,H strongly in L 2 (Q T ). 

After taking a subsequence if necessary, we see that 

E n — > E, a.e. in Qt- 

To show the existence of a weak solution to (1.1)-(1.4), we only need to show 

a(x, |E n |)E n -> |E|)E in ^(Qt). 

As cr(x, s) is continuous on s and E n converges to E almost everywhere in Qt, we 
know 

cr(x, |E n |)E n — > cr(x, |E|)E a.e. in Qt- 

We now show that <j(a;, |E n |)E n is equip-integrable in Qt- We adopt a technique used 
for scalar elliptic and parabolic equations. Let A be any measurable subset of Qt- 
For any large m > 0, 

a(x, \E n )\E n \dxdt 

A 



< a(x,\E n \)\E n \dxdt + / / a(x,\E n \)\E n \dxdt 

J JAf]{\E n \<m} J JAf|{|E„|>m} 

= h + h- 

The assumption on a(x, s) yields 



h<C [1 + \E n \ p \E n \}dxdt, 

■ lAf]{\E n \<m} 
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which can be arbitrarily small if \A\ is small since E n e L p 
On the other hand, 



(Qt). 




IL 



a(x, \E n \)\E n \ 2 dxdt < 



C 



Af|{|E„|>m} 



m 



which is also small if m is sufficiently large. 

This concludes that a(x, |E n |)E n is equip-integrable in Q T . 
It follows by Vitali's theorem that 



a(x, |E n |)E. 



'n 



a(x, |E|)E 



in L\Q T ). 



Finally, we show that (E,H) is a weak solution of (1.1)-(1.4). By multiplying 
Eq.(2.11) and Eq.(2.12) by test functions * and <&, respectively, and then taking 
integration over Q T , after some routine calculations and taking the limit, we see that 
(E,H) is a weak solution to the system (1.1)-(1.4). 

Now we consider the case where a(x, s) is discontinuous on s at some points. 
Without loss of generality, we may assume that a(x, s) has a jump only at one point 
s — 1. In this case a(x, s) is not uniquely defined at s — 1. We shall understand the 
value of a(x, 1) in the following sense: 



where a(x, 1±) represents the right or left limit as s — > 1. 

By the standard approximation, we can construct a smooth approximation se- 
quence <J m (x, s) such that 



Let (E m , H m ) be a solution of (1.1)-(1.4) in which a(x, s) is replaced by a m (x, s). 
By the same energy estimate we see that there exists a measurable function /3(x, t) e 

p+2 

Lp+i (Q t ) such that 



a(x, 1) G [a(x, 1-), a(x, 1+)] 




a m (x, s) is monotonic increasing for all s > 0, 
a m (x, s) = a(x, s), if \s - 1| > ^. 



a m (x, |E m |) -> P(x,t), weakly in Lp+i(Q t ). 
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Define 



A m = {(x,t):l--<\-E\<l + -}, 
m m 



A={(x,t):\E(x,t)\ = l}. 
Since a(x, s) is continuous except at s — 1, we see 

(3(x, t) = a(x, |E|)E, if (x, t) G Q T \A. 

Now it is clear that 

oo 

A = □ An- 

m=l 

Recall that a m (x, s) = cr(x, s) if \s — 1| > ^. It follows that for all (x, t) G A m 

a(x, 1 - — ) < a m (x, |E|) < 1 + — ). 
m m 

Consequently, as m — > 0, 

0"(z, 1-) < /?(x, t) < (7(x, 1+), (x, t) G A. 

Thus, (E,H) is a weak solution of (1.1)-(1.4). 

Finally, we show the uniqueness. Suppose (E,H) and (E*,H*) are two solutions 
of (1.1)-(1.4). Let 

E = E-E*,H = H-H*. 
Similar to the calculation in deriving energy estimates, we find 

sup /[|E| 2 + \H\ 2 }dx+ f f [<t(x, |E|)E — a(x, |E*|)E*] • [E — E*] dxdt 

< o. 

The monotonicity of a(x, s) implies that the second term in the above inequality is 
nonnegative. It follows that 



sup f [|E| 2 + \H\ 2 ]dx < 0. 

Kt<T 



0<t<T- 

Therefore, the uniqueness follows immediately. 

Q.E.D. 
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3. Singular Limit Problem 



In this section we shall show that the solution of (1.1)-(1.4) has a limit as e — > 0, 
which solves Maxwell's equations in quasi-stationary fields, i.e. the system (1.1)-(1.4) 
with e = 0. A weak solution of the quasi-stationary system is defined as in Definition 
2.1 with 6 = 0. 

From now on we denote by (E e ,H e ) the weak solution of the system (1.1)-(1.4). 

Theorem 3.1: The limit of (E e ,H e ) as e — > solves the quasi-stationary system 
(1.1)-(1.4) with e = in the weak sense. Moreover, the weak solution is unique if 
a(x, s) > for all (x, s) G Q x R + . 

Proof: The crucial step in proving the convergence is to show 



where 3(x,t) e Lp+ 1 (Qt). Moreover, as divH £ (x,t) = 0, by the decomposition 
property of H 1 ^) property, after extracting a subsequence if necessary we see that 

H e — > H, strongly in L 2 (Q T ) 



<r(|E e |)E e -> ff(|E|)E, a.e. in Q T as e -> 0. 



The monotonicity of o~(x, s) in s plays a key rule. 

First of all, from Lemma 2.1 and the weak compactness we see 




P+2 



and 



H e — > H, a.e. in Q T . 



Next we show 

J(x, t) — a(x, |E|)E, a.e.inQT- 
We use a monotonicity argument. As a first step, we show 
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Here we adopt an idea from |H|. Let X(t) be a nonnegative smooth function and 

X'(t) < 0,A(0) = 1,A(T) = 0. 
Define an operator L in L p+2 (Qt) 3 as follows: 

L[E] = a(x, |E|)E. 

Since a(x, s) is monotonic increasing in s, then the operator L is monotonic increasing, 
that is, 



<L[E e ]-L[E],E e -E>>0. 



It is clear that 



<L[E 6 ]-L[E],E e -E> 

= < L[E e ],E £ > - < L[E e ],E > - < L[E],E e > + < L[E],E > 



It follows that 



Urn in/ < L[EJ,E E >>< J,E > . (3.1) 

£-+0 

On the other hand, from the system (1.1)-(1.2) we have 

[ T [ a(E £ )\E £ \ 2 X(t)dxdt 
Jo Jn 

= - ( ( [eX{t)E £t ■ E E + XH et ■ HJ dxdt + f f F E £ dxdt 
Jo Jn J Jq t 



lo Jn dt 



|(e|E £ | 2 + |H e | 2 )A(t) 



- m 



Qt 

^(^|E £ | 2 + |H £ | 2 ) 



dxdt + 



F • E £ dxdt 



T 



i|H.|»" 



dxdt H — , 
2Jn 



e|E £ (x,0)| 2 + |H e (x,0)| 2 dx + 



'Qt 



F • E £ dxdt. 



(3.2) 
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Since X'(t) < 0, it follows that 

limsup [ T ( a{\E £ \)\E £ \ 2 X'{t)dxdt 

e^O Jo 

<\ I I \(t)\H\ 2 dxdt + I f \H (x)\ 2 dx + f f E • Fdxdt. (3.3) 
zJo jo 2 Jn J Jq t 

Recall from Definition 2.1 that (E e ,H e ) satisfies the following integral equations: 

[ T f [sE £ ■ ® t + a(x, |E e |)E e • $] dxdt 
Jo Jn 

= [ T [ H £ • (V x &)dxdt + [ [ F-& £ dxdt + [ eE ■ <f>(x, 0)dx, 
Jo Jn J Jq t Jn 

[ T [ [-H e • ^ t + E £ • (V x *)] dxdt = [ \Hq(x) ■ 0)] dx. 
Jo Jn Jn 

Note by Lemma 2.1 that 

i/7 

Jo Jn 



rT 

eE £t ■ <&dxdt\ 

in 



<[ I ! e\E £t \ 2 dxdt] iff e\$\ 2 dxdt 
\Jo Jn I \Jo Jn 



V2 / „ T „ \ 1/2 



^ as e — > 0. 
We take the limit as e — > to obtain 



jf jf [J • $] dzd* = jf jf H • (V x $)d:rdt + | jf F • fcdzd* (3.4) 
/ T f [-H • *t + E • (V x *)] dztft = / \H (x) ■ ¥(x, 0)1 da. (3.5) 

JO Jf7 Jf2 

Now by choosing <3? = A(£)E and * = A(t)H (note that the condition at t = T is 
satisfied since A(T) = 0), we obtain 

' [A(f)J • Eldardi = \ F I [X' (t)\U\ 2 ]dxdt + I [ \H (x)\ 2 }dx +11 E • Fdxd*. 
2 Jo Jn 2 Jn J Jo T 



o 



It follows by (3.3) that 



lim sup < X(t)L[E £ ], E e ><< A(*)J, E > . (3.6) 

e^O 
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Combination of (3.1) and (3.6) yields 

lim < A(t)L[E e ],E c >=< A(f)J, E > . (3.7) 
Consequently, by choosing X(t) properly we have 

lim < L[E e ], E e >=< J, E > . (3.8) 



£^0 



For any vector field W G L 2 (Q T ) f| L P+1 (Q T ), the monotonicity of cr(x,s) in s 
implies 

// [a(x, |E e |)E e - a(x, |W|)W] • [E n - W]dxdt > 0. (3.9) 

We rewrite the above inequality to the following form: 

J J {a(\E £ \)\E £ \ 2 -a(\E £ \)E £ -W}dxdt 

> [ [ [a(x, |W|)W] • [E £ - W]dxdt. (3.10) 
J Jq t 

We take the limit as e — > and use (3.8) for the first term in (3.10) to obtain 

{ J • E - J • W} dxdt 

> f f [u(x, |W|)W] • [E — W]dxdt. 
J Jq t 

Equivalent ly, 

// {[J-<r(a;,|W|)W]-[E-W]}cZa;dt>0. (3.11) 

J J Q T 

Set W = E + 5Y, where 5 > is small parameter and Y G L 2 (Q T ) f| L p+1 (Qt) is 
arbitrary. 

With the above choice of W in the equality (3.10), we obtain 

{Y • [J - a(\E + SY\)[E + 5Y]}dxdt > 0. 

Qt 
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When a(x, s) is continuous in s, then we let 5 — > to obtain 

J(x,t) = cr(x, |E|)E, 

since Y(x,t) is arbitrary in L 2 (Q T ) C\L p+1 (Qt)- 

When a(x, s) has a jump at a point, say, s — 1. Then as in §2 we understand that 
the value of cr(x, s) at s = 1 is 

(t(x, 1—) < cr(x, 1) < a(x, 1+). 

By using the same procedure as in §2, we can derive the above inequality. 

Finally, by taking limit for (2.1)-(2.2) we see that (E, H) is a weak solution of the 
quasi-stationary system. 

To prove the uniqueness, we assume that (E, H) and (E*, H*) are two weak solu- 
tions to the quasi-stationary system. Let E = E — E* and H = H — H*. Then the 
energy estimate implies 

sup / \H\ 2 dx+ f f [a{x, |E|)E-<7(x,|E*|)E*] • \E-¥?]dxdt < 0. 

0<t<TJO. J JQ T 

The monotonicity of a(x, s) in s implies the second term in the above inequality is 
nonnegative. It follows that 

H = 0, a.e.in Qt- 
From the definition of weak solution, we have E = as long as a > 0. 

Q.E.D. 

To conclude this section, we consider two special classes of electric conductivity 
a(x, s) = s p and 

a, if \s\ < 1; 
6, if \s\ > 1, 



a(x, s) 



where < a < b. 

For the first case with F = 0, it is easy to see that 



E = IV x H|"5+TV x H. 
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It follows from Theorem 3.1 that there exists a unique weak solution to the following 
evolution system 

H ( + Vx[|Vx HpifrV x H] = G, (x,t) e Q T , 

subject to the initial-boundary conditions: 

Nx(VxH) = 0, (x,t) 6 dtt x [0,T], 

H(x,0) = H (x), xeQ, 

where G is a known exterior magnetic field. This existence result is not covered in 
[|16|1 . More regularity of the weak solution can be established as in PI . We shall not 
repeat it here. 

For the second case, we see that 

E = -V x H + -F, 

a a 

in the region = {(x, t) : |E| < 1} and 

E = JvxH + V 

b b 

in the region = {(x,t) : |E| > 1}. Note that 

divH(x,t) = divHo(x) = 0, 

we see that 

V x V x H = -AH. 
It follows that H satisfies the parabolic equation: 

H t - - AH = --V x F, (x, t) e Q T 

a a 

and 

H< - j-AH = --V x F, (x, t)eQ$. 
The regularity theory of parabolic equations implies that H(x,t) is smooth in Q^. 
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The interface between Q T and Qj> is defined by 



r = {(:M):|E| = l}, 

which is a free boundary. 

Remark 3.1: We may allow that o~(s) — if |E| < 1 m £/ie above example. In this 
case, one must consider the full system (1.1)- (1.2) in order to define a weak solution. 
However, in this case the uniqueness of the weak solution does not hold. 
Remark 3.2: It is not clear whether or not V is indeed a hyper surf ace in B? x (0, oo] . 
It would be of great interesting to study the smoothness of the interface T and to find 
the free boundary conditions for H. 
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